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"METHOD TO IDENTIFY THE MECHANICAL PROPE RTIES OF A | 
MATERIAL " 

* * * * * 

FIELD OF THE INVENTION 
5 The present invention concerns a method to identify 

the mechanical properties of a material, such as for 
example a bone material, starting from at least a 
digitalized image of a sample of said material. The image 
of the sample to be analyzed is processed to achieve a 

10 numerical model comprising a determinate number of 
elementary geometric elements, hereafter called cells, to 
each of which a parameter to identify the type of material 
is assigned, in order to obtain information on the 
distribution of the mechanical characteristics in the 

15 individual units of the sample, from which it is possible 
to trace both the elastic characteristics and hence the 
resistance, and also the distribution of stresses and 
strain in relation to the application of the loads. 

BACKGROUND OF THE INVENTION 

20 It is known that to diagnose and treat bone 

pathologies, such as for example osteoporosis, it is first 
of all necessary to evaluate the mechanical properties, 
particularly elasticity and resistance to applied loads. 

In order to estimate the resistance and elasticity of 

25 the bone with a good level of accuracy, however, it is not 
enough to know its composition and density. The elastic 
properties of the bone tissue depend to a large extent on 
the architecture of the bone structure in the sample in 
question, as is known, for example, from: Martin, R.B. 

30 (1991) J. Biomech., 24, 79-88, and from: Waud, C.E. et al. 
(1992) Calcif. Tissue Int., 51, 416-418; and mechanical 
resistance is linearly proportionate to the value of Young 
modulus [Hodgkinson, R.J., Currey, D. (1990) P.I.M.E., 204, 
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115-121] . At present, the mechanical properties of the 
bones, and in particular the parameters related to 
elasticity and strength, cannot be determined by mechanical 
tests performed on the patient. Nor, for obvious reasons, 
5 is it appropriate to use bone samples taken from the 
patient on which the tests are to be performed in machines 
suitable for mechanical characterization. On the other 
hand, techniques for processing medical images have become 
very accurate, making it possible to perform digital 

10 modelling of the bone structure on different levels of 
definition, up to the range of some tenths of a micron. 

The numerical methods at present employed for 
modelling and simulating the behavior of the bone structure 
for the purposes of estimating the mechanical properties 

15 have recourse to implementations of the finite elements 
method, known for example from Homminga, J. et al. (2001) 
J. Biomech., 34, 513-517. 

This method, however, has the disadvantage that it 
does not allow to discretize the image easily, that is, to 

20 create a set of geometrical elements able to reproduce the 
structure examined as faithfully as possible and on which 
it is then possible to perform processing and elementary 
quantifications . 

This operation, moreover, is long and laborious, to a 

2 5 large extent entrusted to the ability and experience of an 
operator, since the set of geometric elements has to be 
adapted to the conformation and irregularities of the bone 
structure examined. 

Another application of the finite elements method 

30 provides to make a geometric element correspond to every 
voxel, or elementary volumetric unit, of the image 
acquired; however, this solution has the disadvantage that 
there is a very high number of elements present, such that 
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processing the numeric model generated requires the use of 
computer networks which are not available in normal centers 
of analysis. 

If the resolution, that is, the degree of detail, of 
5 the image is diminished in order to reduce the number of 
elements, the validity of the solution obtained inevitably 
declines . 

In the case of identification of the bone micro- 
structure, another application of the finite elements 
10 method provides to use, as finite elements, rods of 
different length and section, constrained together so as to 
reproduce the trabecular structure. 

This type of modeling generates, on the one hand, a 
number of unknowns low enough to be processed without 
15 particular strategies, but, on the other hand, the result 
obtained is not accurate enough. 

The present Applicant has devised and embodied this 
invention to overcome these shortcomings of the state of 
the art, and to obtain other advantages. 
20 SUMMARY OF THE INVENTION 

The present invention is set forth and characterized 
essentially in the main claim, while the dependent claims 
describe other innovative characteristics of the invention. 
One purpose of the invention is to perfect a method by 
25 means of which it is possible to identify, in an accurate 
manner, the specific mechanical properties of a material, 
such as for example a bone material, in relation to the 
particular conformation of the structure which the material 
assumes . 

30 Another purpose is to perfect a method to calculate 

the mechanical properties of a material irrespective of the 
presence or absence of other heterogeneous materials 
incorporated in the base material. 
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In accordance with these purposes, a method to 
identify the mechanical properties of a material according 
to the present invention comprises a first step to acquire 
at least one digital image of a sample of the material to 
5 be examined, a second step which, for every elementary unit 
of said digital image, acquires at least one predetermined 
parameter indicative at least of the presence or absence, 
and/or the type of matter in the elementary unit itself, a 
third step to construct a grid of elementary geometric 

10 elements or cells able to be associated with the digital 
image, and a fourth step to assign, for each cell, an index 
to represent the characteristics of the type of matter of 
the sample of material which occupies the position of a 
relative cell, in relation to the parameters possessed by 

15 the elementary image units associated with said cell. A 
fifth step defines the mechanical characteristics of each 
cell according to the corresponding index and the 
mechanical properties of the type of matter, and a sixth 
step determines the mechanical properties of said material 

20 according to numerical simulations performed with the known 
Cell Method (Tonti, E. CMES, 2, 237-258) . To be more exact, 
in the course of said numerical simulations, the strains, 
the stresses and/or the distribution of the stresses 
suffered by the sample examined are analyzed, according to 

25 the mechanical characteristics of the cells, defined during 
the aforementioned fifth step. 

In the method according to the present invention, an 
image or a series of images, such as for example a micro-CT 
(computerized micro-tomography) , a micro-MRI (micro- 

30 magnetic resonance), a CAT (computerized axial tomography), 
an NMR (nuclear magnetic resonance) , a digital or 
digitalized radiography or similar, of the portion of 
material examined, is divided into elementary units, voxels 
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in the case of three-dimensional images, or pixels in the 
case of two-dimensional images. 

According to a variant, before being further 
processed, the image is subjected to at least one filtering 
5 process in order to eliminate the distortions due to the 
measuring instrument. Subsequently, for every elementary 
image unit, the possible presence of one or more than one 
of the materials searched for is identified. In other 
words, by means of appropriate calculation mechanisms that 

10 use, for example, threshold algorithms or other similar or 
comparable method, the presence or absence is 
discriminated, in the specific elementary image unit, of 
the material of the sample examined. The data concerning 
the distribution of the mass in the sample is then 

15 processed so as to form an image that can be black and 
white and comprise a suitable number of intermediate shades 
of gray. Each elementary image unit is identified, in this 
case, not only by the three spatial coordinates, but also 
by the predetermined parameter comprising, in this case, 

20 the relative value of the shade of gray. Examining the 
shade of gray in a bone material also allows to 
discriminate the presence of an implant or a prosthesis. 

According to one solution of the invention, above the 
image a grid is defined of some thousands points, or nodal 

25 points, which serve to construct the numerical model for 
the structure calculation forming the cells on which 
processing is performed. The grid of nodes covers the whole 
image, irrespective of the presence or absence, at some 
points thereof, of matter. The distance between the nodal 

30 points constitutes the lower limit of the size of the 
calculation structure under examination, therefore the 
definition of a desired sensitivity imposes a minimum 
distance between the nodes to be used in the specific 
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application . 

Using the base mechanical characteristics, or 
constitutive law, of the material examined as a calculation 
parameter, such as for example elastic-linear, elastic- 
5 plastic, with or without work hardening, isotropy or 
anisotropy or similar, the mechanical characteristics of 
each cell are found also according to the value assumed by 
the corresponding index obtained from the indexes of the 
corresponding elementary image units. 
10 The calculation method thus obtained allows to perform 

an analysis under load by means of the aforementioned Cell 
Method. 

According to a first operative solution, the aforesaid 
Cell Method is applied in the static formulation, although 

15 other formulations could also be applied. 

The Cell Method provides to simulate the application 
of loads on the structure. Moreover, load tests are 
simulated, in the sense of application of forces, pressures 
or torques, or displacements imparted at points of the 

20 structure; the load may be concentrated or distributed at 
points, lines or surfaces of the cells, or on their 
relative combinations, in predetermined directions, in 
order to obtain a plurality of specific displacements of 
the nodal points of the model following the application of 

25 the simulated loads. The static analysis of the 
displacements obtained allows to determine a plurality of 
mechanical properties of the material examined, such as the 
elastic properties of the structure, the distribution of 
stresses and strains, and the zones of maximum stress 

30 inside the material itself, and also the equivalent 
stresses according to different failure theories. 

The method according to the present invention has the 
advantage, compared to the traditional finite elements 
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method, that it automates the division into discreet 
geometric elements, such as cells, without needing to have 
recourse to an excessive number of unknowns and also when 
there is a heterogeneity of materials or when there is no 
5 matter in some points. 

When a grid of nodal points is used, an index is 
attributed to each of them, proportionate to the value of 
the shade of gray of the elementary unit (voxel, pixel...) 
of the digital image. The nodes are then connected, to 

10 create a geometric element or cell, triangular in shape or 
any other geometric shape. A further index is attributed to 
each cell, which is a function of that of its nodes, which 
is representative of the structure of the material that, in 
the image, occupies the position of the cell in the 

15 calculation model. 

In a first embodiment, the cells are determined once 
and for all and entirely cover the whole region, plane or 
spatial, which is examined. 

In a variant, a set of cells is created locally from 

20 each node. One node is considered at a time and the node 
examined is taken as the vertex, while a certain number of 
nearby nodes are used to construct a local complex of 
cells, which all have a vertex in the same node examined. 
This process is repeated for all the nodes. Since a limited 

25 number of cells are considered at a time, the process is 
simpler, especially in the case of a three-dimensional 
model . 

The method according to the present invention can 
easily be implemented in any processor, whether it be an 
30 electronic processor or the processor inside a machine, for 
example a machine for digital radiography. The processor 
can be used directly by a medical operator for the accurate 
diagnosis of the bone properties, or to evaluate the 
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validity of a prosthesis implanted in a patient, for 
example, by means of determining the flow lines of a load 
or quantifying the state of stress and/or strain. 

BRIEF DESCRIPTION OF THE DRAWINGS 
5 These and other characteristics of the present 

invention wi 11 become apparent from the following 
description of a preferential form of embodiment, given as 
a non-restrictive example, with reference to the attached 
drawings wherein: 
10 - fig. 1 is a digital image obtained with a MICRO-CT 
instrument; 

- fig. 2 is a reconstructed image in three dimensions 
starting from a series of digital images like those in 
fig. 1; 

15 - fig. 3 is a processed image of some sections of the image 
in fig. 2; 

- fig. 4 is a complex of cells of a triangular type; 

- fig. 5 is a discreet structure image corresponding to the 
image in fig. 3; 

20 - fig. 6 is a processed image of other sections of the 
image in fig. 2; 

- fig. 7 is a discreet structure image corresponding to the 
image in fig. 6; 

- fig. 8 is a first digital radiographic image; 

25 - fig. 9 is a processed image of a zone of the image in 
fig. 8; 

- fig. 10 is a discreet structure image corresponding to 
the image in fig. 9; 

- fig. 11 is a second digital radiographic image; 

30 - fig. 12 is a processed image of a zone of the image in 
fig. 11; 

- fig. 13 is a discreet structure image corresponding to 
the image in fig. 12. 
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DETAILED DESCRIPTION OF SOME 
PREFERENTIAL FORMS OF EMBODIMENT 
In accordance with a first form of embodiment, with 
reference to figs. 1 to 7, the method according to the 
5 present invention is used to identify the specific 
mechanical properties of a bone sample of about 8 mm in 
diameter and 10 mm in height, taken from a pig's femur and 
subjected to a chemical stabilization process. A digital 
micro-CT image (fig. 1) of the bone portion was obtained, 

10 in synchrotron light, at the SYRMEP line (ELETTRA - 
Trieste) of a 3.54 x 6.72 mm portion inside the sample. 
Using a portion inside the sample allows to avoid improper 
alterations on the edges following the cutting procedure 
with which the sample itself is obtained. The method 

15 provides to assume, as an input parameter, the mechanical 
characteristics of the base bone structure, and 
particularly the constitutive law. The image is then first 
processed by means of appropriate algorithms in order to 
eliminate artifacts generated by the measuring instrument. 

20 The images thus processed comprise hollow regions and full 
regions, divided into elementary units of a three- 
dimensional image, such as voxels, each identified by three 
coordinates and a fourth datum related to the shade. The 
image processed in three dimensions (fig. 2) was sectioned 

25 on four parallel adjacent planes, corresponding in all to a 
thickness of 56 micron, in order to obtain four images of 
the bone sample. To each of these images a threshold filter 
was applied, using the convention of assigning the value 
"0" to the hollow regions and "1" to the regions with bone. 

30 The four images were added, so that each point of the 
image, according to the quantity 0 or 1, assumed a value of 
a shade of gray between 0 and 4. The result can be 
represented graphically with a two-dimensional image of the 
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type shown in figs. 3 and 6. 

On the resultant image a grid of 27 69 nodes 11 was 
deposited (fig. 4} . The average distance between the nodes 
11, depending on the sensitivity of the results to be 
5 obtained, was fixed at 106 micron. The nodes 11 were 
connected to each other so as to form a complex of 5340 
triangular cells 12. An index value to be attributed to 
each cell 12 was then obtained by adding together the 
values of the indexes of the corresponding nodes 11 which 

10 constitute the vertexes thereof and normalizing the result 
to 1. The physical significance of the index of the cell 12 
as thus defined is that it gives an indication of the 
percentage of bone in the cell 12. Fig. 5 shows the result 
of this operation. 

15 The base mechanical characteristics were assumed for 

the bone, such as an elastic-linear isotropic constitutive 
law with a Poisson ratio of 0.3 for all the cells 12 that 
have an index other than 0. The cells 12a with an index of 
0 represent cavities, that is, boneless portions, and do 

20 not possess any mechanical characteristics. In the cells 
12b with an index of 1, a Young module (elastic module) was 
assumed, equal to 1 GPa . In the cells 12c that have 
intermediate index values, the Young module was related to 
the value of the index. For example, if in a cell the index 

25 value is 0.5, the Young module of the cell is 0.5 GPa, and 
so on. Fig. 5 is therefore representative of the values 
assumed by the Young module of each cell 12. Moreover, a 
suitable factor was used to define the porosity of the bone 
structure examined. Given that the index (I ce u) of each 

30 cell 12 represents the quantity of bone material present in 
the cell, the sum of all the indexes divided by the number 
of cells 12 (N cell8 ) indicates the percentage of bone 
structure in the sample. The complement of 1 is therefore 
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For the bone structure examined the porosity factor 
was equal to 60.9%. 

The model of cells 12a, 12b, 12c was then subjected to 
5 a compression test simulated along the long side of the 
image (x axis) . However, it comes within the field of the 
present invention that compression and shear tests can be 
simulated in different directions, or simulations of 
different types of test, for example with one or more loads 

10 concentrated or distributed (pressures) at points of the 
model, also as chosen by the operator, and different 
constraint conditions can be recreated, so as to obtain the 
values of the displacements of the nodes 11 of the model, 
following the application of the simulated loads. In this 

15 specific case, displacements of the nodes on one end (the 
vertical side of the image) of the model of cells 12 were 
constrained in the horizontal direction (x axis) in the 
nodes 11 and in both horizontal and vertical (x and y axes) 
in the lower corner node; a displacement along the x axis 

20 was applied to the nodes 11 of the other end (the vertical 
side of the image) , which simulates the displacement of the 
transverse bar in a test machine. The displacements in 
direction y were left free. 

From the analysis, following the Cell Method, of the 

25 values of displacement, the first datum obtainable from the 
simulated mechanical test was the compression value of the 
apparent Young module along the x axis of the model of 
cells 12 or the bone structure examined, which was 143 MPa. 
The tensions acting in each cell 12 as a consequence 

30 of the load applied in the simulated test were also 
evaluated. To be more exact, the principal stresses were 
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identified and the most stressed points were evaluated with 
respect to a nominal stress value of a structure with zero 
porosity, that is, with all the cells 12 with a unitary 
index. This procedure also allows, for a more accurate 
5 prediction of the risk of fracture, to calculate the 
stresses and identify the zones of greatest risk according 
to two known failure theories, such as a maximum normal 
tension, or Rankin's theory, and a theory of distortion 
energy, or von Mises' theory. 

10 A second portion of the sample was then examined, 

obtained from the same three-dimensional reconstruction of 
the sample (fig. 2), but this time sectioned in a different 
direction (fig. 6) . Given that the mechanical properties of 
the bone structure are very anisotropic, that is, they 

15 depend on the direction in which the load is applied, the 
comparison with the previously examined structure allowed 
obtaining indications on the validity of the system 
proposed. 

The grid of nodes was the same as the one used before, 
20 and also the complex of cells. The discretized structure is 
shown in fig. 7, and in this case too a porosity factor of 
60.9 was obtained. 

A simulation was made for this second structure too, 
performed with the same methodology as the previous one and 
2 5 analyzed with the Cell Method. 

The apparent Young module of the model of the second 
bone structure was equal to 37 MPa, a much lower value than 
that found in the first simulation. 

This result is due to the fact that the bone structure 
30 is arranged in a very different way in the two cases, and 
therefore, even when there is an equal porosity value, 
different results are obtained from the structural 
analysis. Moreover, the trabecula of the first structure 
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are mainly oriented along the axis of compression in the 
simulated test (long side of the image), whereas this is 
not so for the second structure, which in fact has a great 
decline in the mechanical properties. 
5 In accordance with a second form of embodiment, with 

reference to figs. 8 to 13, the method according to the 
present invention is used to identify the specific 
mechanical properties of the bone structure in humans. A 
digital radiographic image (fig. 8) was obtained of a human 

10 hand belonging to a clinically negative subject (no 
osteoporosis) by means of a digital radiography with a flat 
panel produced by General Electric and taken in a hospital 
department. As in the previous embodiment, the method 
provides to assume as an input parameter the mechanical 

15 characteristics of the base bone structure, particularly 
the constitutive law. Inside the radiographic image a zone 
is identified wherein the investigation is performed, in 
this specific case, in the proximal zone of the first 
phalange of the ring finger, although, obviously, many 

20 other parts of the body can be used. 

The image of the zone identified is oriented so as to 
have the horizontal direction, axis x, orthogonal to the 
direction of the finger and the vertical direction, axis y, 
along the phalange. The bone matrix is extracted by means 

25 of applying particular non-linear filters (in this specific 
case a sub- threshold erosion filter, but other 
implementations are possible, such as for example part 
threshold filters, discrimination filters, filters in the 
space and frequencies domain, etc.). After processing, each 

30 pixel of the image in the selected zone has a value of the 
shade between 0 and 1. The result can be represented 
graphically with an image of the type shown in fig. 9. On 
the resultant image a grid of 5662 nodes was deposited. The 



WO 03/082118 PCT/IB03/01127 

» .j 

- 14 - 

nodes were connected together so as to form a complex of 
11050 triangular cells. An index value to be attributed to 
each cell was then obtained by adding together the values 
of the indexes of the corresponding nodes which constitute 
5 the vertexes thereof, the values of the indexes of the 
corresponding baricenters of each cell and the values of 
the indexes of the corresponding middle points of the sides 
of each cell, and by normalizing the result to 1 (but other 
implementations are possible) . The physical significance of 
10 the index of the cell as thus defined is that it represents 
an indication of the percentage of bone in the cell. Fig. 
10 shows the result of this operation. 

The base mechanical characteristics were assumed for the 
bone, such as an elastic-linear isotropic constitutive law 

15 (but other constitutive laws are possible) with a Poisson 
ratio of 0.3 for all the cells that have an index other 
than 0. The cells with an index of 0 do not possess any 
mechanical characteristics. In the cells with an index of 
1, a Young module (elastic module) was assumed, equal to 1 

20 GPa. In the cells that have intermediate index values, the 
Young module was related to the value of the index. For 
example, if in a cell the index value is 0.5, the Young 
module of the cell is 0.5 GPa, and so on. Fig. 10 is 
therefore representative of the values assumed by the Young 

25 module of each cell. 

The cells model was then subjected to two simulated 
compression tests, one along each side of the image (x and 
y axes) . However, it comes within the field of the present 
invention that compression and shear tests can be simulated 

30 in different directions, or simulations of other types of 
tests, for example with one or more loads concentrated or 
distributed at points, lines and/or surfaces of the model, 
also as chosen by the operator, and different constraint 
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conditions can be recreated, and/or the values of the 
displacements of the nodes of the model can be set. 

In this specific case, displacements of the nodes on 
one end (the vertical side y of the image) of the cells 
5 model were constrained in the horizontal direction (x axis) 
in the nodes and in both horizontal and vertical (x and y 
axes) in the lower corner node; a displacement along the x 
axis was applied to the nodes of the other end (the 
vertical side y of the image) , which simulates the 

10 displacement of the transverse bar in a test machine. The 
displacements in direction y were left free. The method was 
then repeated in an analogous manner, constraining the 
horizontal sides and imparting a displacement along the y 
axis of the image. 

15 From the analysis, following the Cell Method, of the 

values of displacement, the datum obtainable from the first 
simulated mechanical test was the compression value of the 
apparent Young module along the x axis of the cells model 
of the bone structure examined, which was 279 MPa. Along 

20 the y axis, the value of the apparent Young module of the 
cells model of the bone structure examined in the second 
simulated test (along y axis) was 322 MPa, confirming the 
method's ability to identify the anisotropy of the bone 
structure examined . 

25 A second digital radiographic image (fig. 11) was then 

examined, of a human hand, belonging to a clinically 
positive subject (with recognized osteoporosis) . 

Proceeding in the same way as previously described, 
the zone wherein to carry out the investigation is again 

30 identified in the proximal zone of the first phalange of 
the ring finger (fig. 12) . The image is oriented so as to 
have the horizontal direction, x axis, orthogonal to the 
direction of the finger and the vertical direction, y axis, 
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along the phalange. 

The grid of nodes was the same as the one used before, 

and also the complex of cells, the filter method and the 

method to attribute the index to the cells - that is, the 
5 whole method was repeated on the new image. The discretized 

structure is shown in fig. 13. 

Two simulations were made for this second structure 

too, performed with the same methodology as the previous 

one and analyzed with the Cell Method. 
10 The apparent Young module of the model of the second 

structure examined in directions x and y of the image was 

respectively equal to 187 MPa and 258 MPa, lower values 

than those found in the first set of simulations, which 

result is due to the different resistance of the two 
15 structures examined, depending on the presence, in the 

second case, of a pathology which entails a great decline 

of the mechanical properties. 

It is clear that modifications or addition of steps 

may be made to the method as described heretofore, without 
20 departing from the field and scope of the present 

invention . 

It is also clear that, although the present invention 
has been described with reference to some specific 
examples, a person of skill in the field shall certainly be 
25 able to achieve many other equivalent forms of the method 
to identify the mechanical properties of a material, all of 
which shall come within the field and scope of the present 
invention. 



